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ABSTRACT 

An important dynamic parameter of barred galaxies is the bar pattern speed, Hp. Among 
several methods that are used for the determination of Hp the Tremaine- Weinberg method has 
the advantage of model independency and accuracy. In this work we apply the method to a 
simulated bar including gas dynamics and study the effect of 2D spectroscopy data quality on 
robustness of the method. We added a white noise and a Gaussian random field to the data and 
measured the corresponding errors in Hp. We found that a signal to noise ratio in surface density 
~ 5 introduces errors of ~ 20% for the Gaussian noise, while for the white noise the corresponding 
errors reach ~ 50%. At the same time the velocity field is less sensitive to contamination. On 
the basis of the performed study we applied the method to the NGC 3367 spiral galaxy using Ha 
Fabry-Perot interferometry data. We found Hp — 43 ± 6 km s _1 kpc -1 for this galaxy. 



Subject headings: galaxies: individual (NGC 3367) 
n-body simulations — techniques: interferometric 

1. Introduction 

The evolution of disk galaxies is strongly influ- 
enced by the main structures within their disks, 
notably spirals and bars. Bars are a fundamental 
component of mass distribution in spiral galaxies 
because they may contain a large fraction of the 
disk. Bar kinematics is different from the one of 
the axisymmetric disk, which makes the study of 
bar parameters important, especially the bar pat- 
tern speed (ftp). Bar properties are tightly related 
to the mass distribution of the host galaxy. The 
shape and structure of the bar, such as vertical 
bending (buckling ), is believed to dep end on the 
stage of evolution (| Athanassoulall2003l ) . Moreover, 
Hp is one of the fundament al param e ters in driv- 



galaxies: kinematics and dynamics — methods: 

cation of theoretically predicted resonances (Lind- 
blad resonances, corotation), using the rotation 
curv e to extract the periodic mot i on of stars and 
gas (|Tullvl Il974l ). IContopoulosl (|l98fj ) demon- 
strated with orbit calculations that the corota- 
tion radius should be located at the end of a self- 
consistent bar. However, it was shown from a sur- 
vey of early-type galaxies, that indeed the bar ends 
rather bet ween the inner 4 : 1 res onance and the 
corotation (jElmegreen et al. I ll996h . 



An alternative method consists of matching nu- 
meric al models to the observed data (jRautiainen et al 
20051 ). Two other method s are based on the analy 



ing the evolution of bars. iNoguchil (|1987I ) finds 
from numerical simulations that tidally-induced 
bars rotate slowly in comparison with sponta- 
neous bars. Thus, bar pattern speed may help to 
discriminate between spontaneous bars and tidal 
bars. 

One of the first methods applied to bar pat- 
tern speed determination is based on the identifi- 



sis of spiral density waves. lElmegreen et al. (1992) 
apply a computer algorithm to extract various 
types of symm etries from galaxy images, while 



Canzianl (|l993l ) points out the difference in the 



global appearance of the residual velocity field of 
a spiral galaxy inside and outs i de th e corotation 
radius. iTremaine fc Weinberg! ( 1984 ) formulated 
a method to measure the bar pattern speed that 
is independent from any spiral density wave the- 
ory. The principal assumption of the method is 
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that the surface density of the tracer of the grav- 
itational potential (e.g., old stars) satisfies the 
continuity equation, i.e., there is no significant de- 
struction or creation of the tracer over a dynam- 
ical time. This method allows the measurement 
of the bar pattern speed with two observational 
quantities: the surface brightness of the tracer 
and the velocity of the tracer along the line of 
sight. So far the method wa s succ e ssfully ap- 
plied to some twenty galax ies. iKent (1987) and 
Merrifield fe Kuiiken ( 1995f ) apply this method to 
the st ellar bar comp o nent o f the SBO galaxy NGC 



936. Gersscn et al 



(1999) apply the Tremaine- 
Weinberg (TW) method to the stellar component 
of the galaxy NGC 4596. The same authors use 
the method for the galaxies NGC 271, NGC1358 



ESQ 281-31 and NGC 3992 dGerssen et al.l ll999). 



Aguerri et al. ( 20031 ) applied the method for the 
galaxies ESQ 139-G009, IC 874, NGC 1308, NGC 



1440 and NGC 3412. iDebattista fc Williams! 



()200lh use the TW method for NG C 7079 and 



for NGC 1023 (De battista et all 2002) 



Other authors derived the bar pattern speed 
applying t h e TW method to a gaseous tracer. 
Westpfahll (|1998T ) finds the bar pattern speed 
for M81 empl oying HI as a tracer . Using CO 



observations, IZimmer et al.l (|2004l ) employ the 



method for the g alaxie s M51, M83 and NGC 6946. 
Rand fc Wallinl (|2004h consider the galaxies NGC 
1068, NGC 3627, NGC 4321 (M100), NGC 4414, 
NGC 4736 and NGC 4826, and measure the bar 
pattern speed with CO as a tracer. The applica- 
tion of this method to the gaseous phase is more 
delicate because of the assumptions of the method 
itself. In general, the ionized gas, the CO and HI, 
do not satisfy the continuity equation over an or- 
bital period and do not trace the gravitational 
potential. However, assuming that the gas is con- 
tinuously distributed along the pattern, and that 
the luminosity- weighted mean is a valuable indica- 
tor of the mean mass distribution and, neglecting 
the internal kinematics of H II regions, some au- 
thors successfully applie d the TW method to the 
ionized gaseous phase. Hernandez et al. ( 20041 ) 



employ ionized hydrogen 2D velocity fields of four 
barred galaxies NGC 4321, NGC 3359, NGC 6946, 
NGC 2903 and also of M51 to measure the pattern 
speed. The same authors build numerical simu- 
lations to investigate the possibility of using the 
gaseous component and present an application 



for NGC 4321 galaxy (jHernandez et all l2005fh 
(2006j) have successfully applied 



Emscllcm ct al 



TW method for bar pattern speed determination 
of NGC 1068 using Fabry-Perot Ha map, and 
also reproduce with Nbody+S PH models a num- 
ber of its observed properties. iFathi et alj (|2007l ) 
also with two dimensional ionized hydrogen kinc- 
matic s determine the bar p attern speed of NGC 
6946. Beckman et al. (2003) apply the method to 
a sample of nine galaxies (NGC 3049, NGC 4294, 
NGC4519, NGC 5371, NGC 5921, NGC 5964, 
NGC 6946, NGC 7479, NGC 7741) using ionized 
hydrogen 2D kinematic data. 

In this article we employ the TW method to 
study the errors in the parameters that may affect 
the pattern speed determination. We test to which 
parameters the method is sensitive with numerical 
simulations and apply this method to the barred 
galaxy NGC 3367. 

2. Numerical test of the TW method 



The Tremaine Weinberg method uses two ob- 
servable quantities: the luminosity weighted ve- 
locity and luminosity weighted density determined 
along a thin stri p (aperture) parallel to the m ajor 
axis of the disk (jTremaine fc Weinberg Il984h . If 
a galaxy is centered at the cartesian coordinates 
such that the major and minor axes are aligned 
parallel to the x and y axes, respectively, then the 
ratio of intensity-weighted velocity and intensity- 
weighted position gives the angular velocity: 



Q p sin i = 



(V(x)) - V 6 



«!/« 



(x) - Xq 



(1) 



Here V sys - is the systemic velocity, and xo is the 
position of the kinematic center of the disk along 
x axis. The ftp is also corrected due to inclina- 
tion i (i = corresponds to face-on disk) . This 
formulation of Merrifield fc Kuijkenl ( 19951 ) allows 
a more accurate evaluation of ttp because the er- 
rors in the dynamical center and systemic velocity 
determination are reduced. Thus, estimating the 
above expression for several apertures, one may 
plot (V) vs. (x) and obtain the averaged value of 

n P . 

One of the conditions required by the TW 
method is that the tracer should satisfy the con- 
tinuity equation. The old stars in SBO galaxies, 
for example, survive long enough to trace the spi- 
ral pattern potential. The main concern regarding 
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the validity of application of the TW method to 
gaseous tracer is that the ISM, composed mostly 
of molecular and atomic hydrogen, is not able to 
trace the potential for a long period of time due 
to short timescale processes such as phase tran- 
sitions, cooling, etc. Nevertheless, it was success- 



fully applied to CO emiss ion lines (jRand fe Wallin 
2004 : IZimmer et "ah 20041) and to the Ha emission 
line (IHernandez et al.ll2005t lEmsellem et alJ 2006; 
Fathi et al.l l2007t) . As was shown by IWestpfahl 
(|l998l ) addition or substraction of the tracer ma- 
terial has no effect on calculation of the instanta- 
neous pattern speed, and thus, it is expected that 
processes o f star formation and fee dback may be 
neglected. iHernandez et al.l (|2005l ) have demon- 
strated by hydrodynamical simulations that the 
application of the method to ionized gaseous phase 
is reliable under some assumptions. The problem 
arises rather in the quality of the image and the 
velocity field as affected by regions of star forma- 
tions, dust obscuration, and local gas motions. An 
argument in favor of applicability of the method 
is that the Ha image also contains the stellar con- 
tinuum emission. This should alleviate the prob- 
lem of a patchy monochromatic image, making it 
smoother. All these arguments are still have to be 
carefully verified, and here we assume as a work- 
ing hypothesis that application of the TW method 
to the Ha kinematic data is valid. 

As it was shown by Debattistal ( 20031 ) the 
method is very sensitive to errors in the determi- 
nation of the position angle (PA) o f the major axis. 
On th e other hand, as shown by iRand fe Wallin 
(|2004l ). the molecular gaseous component may 
produce a non-zero pattern speed even in the ab- 
sence of any clear wave pattern due to dumpiness. 
The latter point requires that the determination of 
Q, p should be done w i th mu ch caution. In particu- 
lar, IHernandez et al.l (|2005l ) have shown by means 
of hydrodynamical simulations that the regions of 
shocks and the zones outside the bar should be 
avoided. The effects of the inclination, bar ori- 
entation, angular resolution, and the uncertainty 
were already in v estiga ted in numerical models by 
Rand fc Waliir] (|2004l ). 



2.1. Numerical models 

In order to test the influence of errors in inten- 
sity and 2D velocity fields on the determination of 
the pattern speed with the TW method we per- 



formed numerical simulations of a bar formed in 
a hydrodynamic (Nbody+SPH) and collisionlcss 
(Nbody) disk galaxy. 

For this purpose we prepared a bar unstable 
disk galaxy consisting of total 3 x 10 6 particles 
sampling an expon ential stellar and gas disks, and 
a Hernquist halo ( Hernquistl 1990). The verti- 
cal structure of the disk is described by isother- 
mal sheets with a constant scale height of 280 pc 
for stars and 80 pc for gas. Gas dynamics obeys 
the isothermal equation of state with a tempera- 
ture T= 10 4 K. The processes of star formation 
and feedback were not included. The model was 
constructe d using a t echniq ue similar to that de- 
scribed bv lHernquistl (j.1993) and t he numerical pa 



rameters were chosen according to lGabbasov et al 
(2006). The galaxy model parameters are sum- 



marized in the Table [TJ The simulation was per- 
formed with GADGET 2 code on HP CP 4000 
cluster (KanBalam) at DGSCA-UNAM. 

We estimate the bar pattern speed as follows. 
At each snapshot we determine the orientation of 
the principal axes of the inertia tensor of the bar. 
Then, we draw the bar position angle as a func- 
tion of time, and the f2p is obtained by numerical 
differentiation, Clp = d<j)/dt. The bar appears at 
a time t ~ 2 Gyrs and at t — 3 Gyrs reaches its 
maximum length of ~ 7 kpc. For t = 4 Gyrs most 
of the gas is transferred to the center of the disk. 
During this period, the angular velocity of the bar 
remains roughly constant, 41 km s _1 . Further evo- 
lution leads to the depletion of the gas particles 
along the bar, and after about 5 Gyrs, the gas is 
rather located in the center and in the spiral arms. 
For our analysis we took the snapshot at t = 3.75 
Gyrs when the bar is oriented ~ 40° from the ver- 
tical y axis. We added some systemic velocity to 



the velocity field and inclined the disk by 



Table 1: Galaxy model parameters. 



30° 



Parameter 


Gas 


Disk 


Halo 


M (1O 1U M ) 
r (kpc) 
N (10 6 ) 

£ (pc) 


0.4 
3.3 
0.2 

80 


3.96 
3.3 
0.8 
80 


22 
6.6 
2.0 
160 



Note. — Here M is the mass of the component, ro - the scale radius, N - the num 
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such that the disk major axis is aligned with the x 
axis. We compute the projected on the sky plane 
2D surface density and the velocity fields on uni- 
form cartesian grids of 204 x 204, see Fig. [TJ We 
apply the TW method to simulated stellar and 
gas bars and compared O p with the pattern speed 
obtained from the simulation. A difference of less 
than 1.0 km s _1 was found for the stellar bar. The 
gas pattern speed was overestimated by ~ 7 km 
s -1 , because of the surface density weighting er- 
rors. We also applied the method to the galaxy 
at t = 1 Gyrs (early stage of evolution without a 
bar) and obtained a zero slope as expected. 



kpc 



120 




0.0 0.2 0.4 0.6 



Fig. 1. — Projected surface density of the gas, 
with the corresponding velocity contours oveplot- 
ted. The surface density is normalized by the max- 
imum value. The contours are separated by 40 km 
s _1 and dashed lines show some of the slit posi- 
tions (one of ten). 



2.2. PA variation 

First, we tested the sensitivity of the bar pat- 
tern speed to variations of the position of the ma- 
jor axis of the simulated disk galaxy, which could 
be interpreted as errors in the determination of 
PA. Here, the variation in the position angle <5PA, 
is the angle between the disk major axis and the x 
axis of the fixed reference frame. The resulting Op 
is presented in Fig. [5] For the considered range 




SPA (degrees) 

Fig. 2. — Dependence of Hp on the PA variation. 



of SPA and bar orientation the Op changes al- 
most linearly. As the bar tends to become aligned 
with the y axis the pattern speed reduces. From 
this plot one may observe that an uncertainty of 
SPA ±5° gives rise to an error in Op of ±15 km 
s _1 kpc _1 for gas bar. 

As an example, in Fig. [3j we plot the (x), {V}, 
and in Fig. 3] the normalized mean surface den- 
sity (E) /Y>max vs. y for 40 apertures of width 
0.12 kpc about the kinematic center located at 
yo = 2.45 kpc. Also in Fig. |4] shown is the re- 
sulting (V) to (x) ratio. In order to avoid a dis- 
continuity in center of these plots we have ex- 
cluded from the analysis the central four pixels 
with very high gas density. The plots are shown 
for 5PA= 0, -10°, +10° (solid, dotted, and dashed 
line, respectively). 

As seen from Fig. [3J the weighted average ve- 
locity is the most sensitive quantity in the plot, 
and the difference comes mainly from the ends of 
the bar, while at the center the slope varies slowly 
for all three cases. At the same time the slope of 
(x) changes in opposite direction than (V), affect- 
ing strongly the Op determined. The decaying 
intensity and velocity profiles in the first and the 
last kiloparsec are due to inclusion of the apertures 
containing spiral arms at the ends of the bar. This 
is observed as a characteristic Z-shape, or a loop- 
shape, in the (V) vs. (a;) plot in Fig. 0] instead 
of a straight line. The weighted averaged inten- 
sity also changes the position and the maximum 
of its shoulders, being nearer to the center and 
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higher for closer alignment of the bar with the y 
axis. Unfortunately, the restriction in y does not 
reduce significantly the errors produced by varia- 
tion of PA. 

2.3. Disk thickness 

According to the original formulation of the 
method, it is valid only for an infinitely thin disk. 
However, galactic disks generally have a certain 
thickness and line of sight velocities also contain 
the projected component of velocity normal to the 
disk (Z component, v z cosi). This component be- 
comes important as the inclination approaches 
degrees, i.e., the galaxy becomes face on. For a 
disk symmetric in the Z direction the vertical mo- 
tions cancel out each other, and the net contri- 
bution is zero. However, the presence of an asym- 
metry in vertical structure (such as bar buckling) 
will also introduce asymmetry in the velocity field. 
We have verified this argument by including the 
projected Z component of velocity in a simulated 
bar velocity field. For our Nbody+SPH model, 
evolved for 5 Gyrs, we found no significant dif- 
ference in Op for both gas and stellar bars. We 
also prepared and ran a pure collisionlcss model, 
incrementing the stellar disk mass by the mass of 
removed gas disk. We traced the evolution of the 
bar up to 10 Gyrs, and quantify the bar buckling 
as a change of r.m.s. of the vertical velocity dis- 
persion of disk particles, a z . Due to heating of the 
stellar disk the vertical dispersions grow linearly 
until the bar begins to buckle at t m 4.5 Gyrs. Af- 
ter the buckling, the r.m.s. dispersions are rapidly 
increased by ~ 10% and this is reflected in a reduc- 



tion of the pattern speed slowdown rate, as shown 
in Fig. [5J The determination of the pattern speed 
becomes more robust as the spiral arms cease, and 
the curve that does not include the Z component 
is much smoother, except for some spikes. In con- 
trast, inclusion of the Z component clearly affects 
the instantaneous Qp determination due to buck- 
ling, but the errors are quite small (a few km s _1 ). 
Thus, for old stellar bars, where buckling is strong 
and the velocity dispersions are high, the vertical 
motions could make an important contribution to 
the observed velocity field introducing errors in 
the flp determination by the TW method. 

2.4. Data quality 

Next, we checked the robustness of the TW 
method to the introduction of an artificial noise 
to the surface density and velocity fields. Such 
noise is naturally present in the observed data. 
First, we add a white noise background to our 
simulated surface density field to represent the 
emission inhomogeneities. The white noise back- 
ground was created by assigning to each pixel a 
random value within (0, £ pe rt) in case of inten- 
sity field and (— V^> e rt) +^p e rt) m case of velocity 
field. Here, S pcr t and Vp er t are the maximum val- 
ues of perturbation of the surface density and ve- 
locity field, respectively. The results are shown in 
Fig. [6] As we increment the amplitude of the per- 
turbation relative to the corresponding maximum 
value of the surface density (S per t/S ma x) or the 
velocity field (|V^ e rt|/|Knax|), the errors become 
dominant. Indeed, for the ratio E pert /E max > 0.3 
(signal to noise ratio < 5) for white noise, the 
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With Z component 




time (Gyrs) 

Fig. 5. — Pattern speed of a collisioness bar determined with TW method, including and excluding the Z 
component (red thin and black thick curve, respectively). 
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Fig. 6. — Variation of flp with the degree of con- 
tamination of intensity (top) and velocity (bot- 
tom). 



determination of the pattern speed is not possi- 
ble anymore because a clear slope cannot be es- 
tablished. As a next step, we test the sensitiv- 
ity of the method to the Gaussian noise. Besides 
the observational noise it simulates local gas in- 
homogeneities in the disk. We create a random 
Gaussian field with spherical (8 = 10) correla- 
tion function and also change its amplitude rel- 
ative to the maximum value of the signal (see 
Fig. [7]) . For this purpose, we used the software 
developed by iKozintsev fc Kedem ( 2000h . When 
added to the intensity field, the signal is also un- 
derestimated, although not so strongly as in the 
case of white noise. Concerning the velocity field 
contamination, one can note that the determi- 
nation of the flp is possible, even reaching the 
ratio |Vpert|/|Knax| ~ 1- However, in the case 
of the Gaussian random field the errors influence 
stronger the resulting (V) vs. (x) slope for the 
gas bar. Note, that the pattern speed is either 
underestimated or overestimated for different bar 
positions and, in case of the Gaussian noise, for 
different realizations of the field. 

In addition, we investigate the effect of a clip- 
ping procedure (i.e., imposing the inferior limit on 
the surface density) and finding that the effect on 
the resulting pattern speed has rather non-linear 
trend as can be seen in Fig. [8] The surface densi- 
ties of both stellar and gas components were nor- 
malized by their respective maximum values. The 
most sensitive part of the curve is due to clipping 
of up to a few percent of the maximum surface den- 
sity, which corresponds to the density of the faint 
spiral arms and the disk that leads to an overesti- 
mation by more than 50% in the pattern speed (cf. 
surface density distribution in Fig. [T] and Fig. 0|. 
A roughly flat region of the curves where the bar 
size due to successive clipping diminishes and be- 
comes thinner, the flp of both gas and stellar bar 
is reaching 65 km s _1 kpc -1 . A rapid drop in 
the curves corresponds to the clipping when the 
bar is not continuous anymore, but rather looks 
like two bright patches, after which the determi- 
nation of the pattern speed is barely possible. This 
behavior holds for several snapshots we have an- 
alyzed between 2 and 4 Gyrs of evolution. The 
same mechanism should be responsible for a rapid 
change in flp for small S pcrt values added to the 
intensity field. Perturbations of the surface den- 
sity of the the same amplitude as of the faint disk 
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and spiral arms lead to underestimation as shown 
in Fig. [6] Thus, we may conclude that the role 
of the spiral arms in the application of the TW 
method is very important. A s imilar conclusion 
was reached bv lRand &: Wallin (2004) who found 
that too much clipping could violate the principal 
assumption of the TW method and lead to incor- 
rect pattern speed determinations. On the other 
hand, we found that the restriction of the maxi- 
mum surface density (clipping "from above" ) pro- 
duces the bifurcation in the center of (V) vs. (x) 
plot (not shown here), creating an effect of pres- 
ence of two different pattern speeds. The size of 
the bifurcated area increases with the diminishing 
of the upper limit of the surface density. If these 
points are omitted from the linear fit, the resulting 
pattern speed is not affected. 




Fig. 8. — Variation of Dp with the degree of clip- 
ping of the surface density. 




tion lines whose rough emissivity is given by 
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Fig. 7. — Random Gauss field used for data qual- 
ity study. 

The observations measure the intensity of Ha 
emission line rather than the gas surface density. 
The results above were analyzed only for the sur- 
face density field of gas and stars derived directly 
from the simulation, but not for the intensity field 
produced by ionized gas. The original TW method 
is formulated for the stellar component, where it 
is assumed that the disk surface density is directly 
proportional to the number of stars per unit area. 
However, in the case of the gas, the tracer is sen- 
sitive to the square of the density because the re- 
gions of ionized hydrogen will produce rccombina- 



2ji-l/2 



(2) 



If we assume that the gas is isothermal everywhere 
in the disk, then the emission map will be a func- 
tion of gas density squared only. In this case, 
our simulations show that Dp is overestimated by 
roughly 30%. Nonetheless, when the emission map 
is clipped at the bar level (i.e., masking the spiral 
arms and disk) the value of Qp is not changed. 
The same result was found for the the observed 
Ha emission map, presented in the next section, 
if transformed into the surface density of the gas. 
The clipping is important in this case because the 
transformation given by equation [2] mainly affects 
the regions with high density gradients, i.e., tran- 
sitions between bar and disk. 

These results establish the validity margins of 
the TW method for our models. Although they 
were determined for a simple galaxy model, they 
give us the general clues on parameter dependen- 
cies. Complete hydrodynamical simulations, in- 
cluding star formation, supernovae feedback and 
radiative cooling will be presented elsewhere. 

3. Application to NGC 3367 

NGC 3367 is classified as a SBc type barred 
galaxy and it is considered as an isolated galaxy 
at a distance of 43.6 Mpc (Ho = 75 km s _1 
Mpc -1 ), located behind the Leo Spur group 
of galaxies. This galaxy has a remote neigh- 
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bor, NGC 3419, at a distance of 900 ± 100 kpc 



( Garcia-Barreto fc RosaddlioOlh . At a radius of 
10 kpc from the nucleus there is an optical struc- 
ture consisting of several Ha knots that rese m- 
bles a bow-shock (Garcia-Barreto et al.lll996al lbl). 



iGarcfa-Barreto et al.l ( 1998 ) report Very Large 
Array (VLA) observations with an angular resolu- 
tion of 4.5" at 1.46 GHz and find radio continuum 
emission from two lobes that extend up to 6 kpc 
outside the plane of the d isk, and a weaker emis- 
sion from the same disk. Gioia et al. I (|l990h also 
reported soft X-ray emission from this galaxy. 

3.1. Ionized hydrogen kinematics data 
analysis 

We use the ionized hydrogen data cube from 
Fabry-Pe rot interferometry obser v ations pub- 
lished in Garcia-Barreto fc Rosadol ( 2001 ). The 
data cube has dimensions 512 x 512 x 48, a final 
image scale of 0.58" p ixel -1 , and spectral sam- 



pling of 19 km s . iGarcia-Barreto fc Rosadol 



(|200ll ) used an interference filter with a central 
wavelength of 6620A, and a narrow band (30A) m 
order to isolate the redshifted Ha emission of this 
galaxy. The authors made the calibration with a 
neon lamp centered at 6598. 95A. The exposition 
time for each channel was of 120 s. No attempt 
was made by the authors for absolute calibration 
of the Ha emission. 

The data were already reduced with the AD- 
HOC^ 

software, and we use it to obtain the in- 
tensity Ha monochromatic map, the radial veloc- 
ity field, and the continuum map. The Ha im- 
age we use for our further analysis was obtained 
by sum of 20 channels, without substraction of 
the continuum. The velocity field was obtained 
by finding the barycenter of the Ha profile peaks 
for each pixel. The photometric center was deter- 
mined from the continuum map by looking for the 
brightest pixel, whic h in this case coincides with 
the k inematic center (jFuentes-Carrera et al.ll2004 , 
2003). We trimmed these maps to 204 x 204 to ex- 
clude external sources and employ a spectral and 
spatial smoothing with a Lorentzian function with 
a FWHM of 3 pixels (~ 1.7" compared with the 
seeing ~ 1.2"), more adapted to the instrumen- 
tal function of a Fabry-Perot. The final Ha image 



and the radial velocity field are shown in (Fig. [9]). 

We use the two dimensional radial velocity 
field to obtain a rotation curv e for NGC 3367. 
Garcia-Barreto fc Rosadol ( 200ll ) already obtained 
a rotation curve for this galaxy using the AIPS 
package. For this work we wish to explore the 
variation of the parameters of the rotation curve 
in order to optimize them because they are of fun- 
damental importance for determining the bar pat- 
tern velocity. For this purpose we use the AD- 
HOCw software. These parameters are: inclina- 
tion i, position angle PA of the major kinematic 
axis, systemic velocity V sys t and kinematic cen- 
te r. We obtain values similar to th ose published 
in Garcia-Barreto fc Rosadol (|200l[ ) with the dif- 



ference that our systemic velocity and position an- 
gle are slightly higher: V sys t — 3164 ± 10 km s _1 
instead of V syst = 3030 ±8 km s" 1 , PA = 60° ±5° 
instead of PA = 51° ± 3°. Our kinematic center is 
located at R.A.= 10 h 46 m 35 a and Dec.= 13°45'00" 
(J2000.0), which are the same v alues reported by 
Garcia-Barreto fc Rosadol (|200ll ). The inclination 
with r espect to the plane of the s k y is i = 30° ± 5° 
as in IGarcia-Barreto fc Rosadol (l200lh . As re- 
ported bv lGarcia-Barreto et al.l (|20071 ) the stellar 



1 http://www. oamp.fr/adhoc/adhocw. html developed by J. 
Boulcsteix. 



bar length is 32" (6.7 kpc) and is oriented 15° ± 5° 
from the kinematic major axis. 

The rotation curve was determined by averag- 
ing velocities in two sectors along the kinematic 
line of nodes. Once obtaining the full rotation 
curve of NGC 3367, we average the rotation curve 
over both sectors. We fit the mean rotation curve 
using a weighted asymptotic regression model. 
The rotation curve of NGC 3367, the averaged ro- 
tation curve and the fit are shown in Fig. 1101 

4. Application of the TW method to NGC 
3367 

We apply the Tremaine- Weinberg method to 
NGC 3367 in order to measure the bar pattern 
speed. We built an IDL program that calcu- 
lates the intensity-weighted velocity (V) and the 
intensity-weighted position (a;) for each strip along 
the kinematic minor axis of the galaxy. Here the 
Ha intensity serves as a weighting function under 
the assumption that the surface density of the disk 
is proportional to the Ha intensity (Edisfc oc Ijja)- 

The TW method involves many parameters 
among which we consider the position angle of 
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Fig. 9. — Ha image (left) and radial velocity field of NGC 3367 with contours levels separated by 20 km 
s _1 (right). The contours on the Ha image depict the bar region and are shown from 280 to 300 with step 
of 5 in arbitrary units. Also shown on the left panel is the position of the major axis and the photometric 
center. The north direction is on top and the east is on the left. 
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Fig. 10. — Full rotation curve of NGC 3367 (left panel) and averaged rotation curve (right panel, solid line) 
with overplotted fit (dashed line). See the text for further details. 
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the kinematic major axis (PA), the minimum and 
maximum of the surface brightness image, posi- 
tion of the kinematic center in pixels, systemic ve- 
locity, inclination, the range of integration along 
the major axis and the location and number of 
apertures along which wc calculate the quantities 
(V) and (x). In order to explore the errors in the 
fip determination associated with the uncertain- 
ties in the parameters we vary each of them to 
build a range of values. We assume that within 
each determined interval there is a subset of val- 
ues where flp is trustworthy if the plots of the 
weighted mean velocity and position show little 
scatter such that the slope in the (V) vs. (x) plot 
is well fitted and passes though the origin. On the 
other hand, in our case the curves of intensity- 
weighted velocity and position should start from 
negative values and move smoothly towards pos- 
itive ones. By varying the full set of parameters 
we found that four of them are actually relevant. 
These are the PA, the minimum of the H Q im- 
age, the length of the slit along the kinematical 
major axis and the range of variation along the 
minor axis of the galaxy. The errors in V sys ti and 
kinematic center have smaller effect on f 2p, as was 
already noted by iMerrifield fc Kuiiken ( 1995 ). 

Thus, we begin with the variation of PA. We es- 
tablish the origin of the Cartesian coordinate sys- 
tem on the photometric center of the H a image of 
NGC 3367, such that the x axis is aligned with 
the disk kinematic major axis. First, we vary the 
PA within the range of errors determined from the 
ADHOCw package (±5°) in both north-east di- 
rection (positive) and in the south-east direction 
(negative). From Fig. [Tl] it can be noted that a 
variation of <5PA by ±5° leads to errors in ftp of 
±6 km s -1 . We extend this range up to 19° in 
order to compare it with the behavior of Qp ob- 
served in the simulation. The result presented in 
Fig. Qj] shows that for i5PA> the curve is al- 
most linear and roughly similar to the Fig. [2l As 
the bar becomes aligned with the kinematic major 
axis (<5PA< 0), the Op approaches zero faster. A 
similar conduct is found for the simulated bar ori- 
ented in the same way as the observed one (~ 15° 
from the major axis). 

In order to select only those emission regions 
near the bar, we mask the H a image by chang- 
ing the minimum of the image. Clipping to a 
given minimum H Q intensity is necessary because 




-10 o 

5PA (degrees) 

Fig. 11. — Variation of ftp with position angle of 
the disk of NGC 3367. 



if the entire image is taken the plot of (V) vs. 
(x) is too noisy and it is difficult to establish 
a clear fit. As shown in the previous section a 
sufficiently noisy background could strongly af- 
fect the results. The masking procedure helps 
to hide bright knots within the interarm regions 
and allows to achieve higher sign al to noise ratio. 
Addit ionally, as demonstrated by iRand fe Wallin 
([20041 ) . intensity clipping helps to improve the re- 
sults by removing the scatter produced by clumps, 
but at the same time it could also remove part 
of the bar pattern. Fig. [12] shows the variation 
of Qp with the minimum of the H Q image. The 
range of varied intensity corresponds to emission 
that roughly traces disk, spiral arms, and the bar 
region, except the bright bulge. When compared 
with the Fig. [8] a similar behavior of the pattern 
speed due to clipping can be observed. Within 
the bar region (excluding disk and arms) the er- 
rors due to clipping are ±5%. The small varia- 
tion of ftp within the clipping range 210 — 240 is 
probably due to the dumpiness of the H Q image 
around the bar. We also examined the clipping of 
the maximum intensity and found that it did not 
significantly alter the results. 

Strong emission along the spiral arms represent 
a great source of noise in bar pattern speed de- 
termination that one needs to isolate in order to 
correctly perform the weighting procedure. The 
emission intensity of spiral arms is much stronger 
than that of the bar, and (x) and (V) will be bi- 
ased by the arms. For this reason we resort to con- 
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Fig. 12. — Variation of £lp with the minimum of 
intensity of the H Q image. 

straining the integration area uniquely to the bar 
zone. We restrict the range of integration along 
the kinematical major axis to avoid the region of 
the eastern arm. This is an extended and promi- 
nent zone that biases significantly the signal to 
noise ratio of the intensity near the north-east side 
of the bar (Fig. [9]). Since we are interested in the 
bar pattern speed, avoiding spiral arms that may 
have a pattern speed distinct from that of the bar 
would also improve the signal. The fourth impor- 
tant parameter is the range along the minor axis 
of the galaxy, i.e., the number of apertures. If 
this range is not limited, the bright spiral arms of 
the galaxy, in particular the one located at north- 
west, enter in the computation and the bar pattern 
speed cannot be determined. 

After analyzing the parameter variation we de- 
termine the bar pattern speed of NGC 3367. The 
parameters mentioned above were constrained in 
the following sequence: limiting the integration 
area, limiting the intensity, and varying the PA. 
The final masked H a image and the region where 
the flp was determined is shown in Fig. 1131 The 
averaging of quantities in equation (JT|) is done 
along strips of one pixel width, totalling 18 slits 
parallel to the kinematic major axis within the 
box. 

We find for NGC 3367 a value of fl P = 43±6 km 
s _1 kpc -1 . This value is consistent with the v alue 
reported in Garcia-Barreto &: RosadJ ( 2001 ). In 
Fig. [14] we show the mean weighted position and 
the mean weighted velocity as a function of the 




Fig. 13. — Masked H Q image and the area of in- 
tegration used for the final Hp calculation. Some 
of the apertures used in the TW method were also 
marked. Box units are in arcseconds, and inten- 
sity is in arbitrary units. 

position along the minor axis of the galaxy. Both 
curves show an almost linear trend, except that 
(x) shows flattening after y — 11 pixels (1.35 kpc). 
This is probably due to a characteristic shape of 
the bar, both parts of which are not symmetric 
and which appears to be "broken" . This asymme- 
try is also observed in broadban d ultraviolet im- 
ages ( Garcia-Barreto et al.l I2007T ) . In Fig. [15] we 
show the mean weighted velocity as a function of 
the mean weighted position and the mean inten- 
sity as a function of the position along the minor 
axis. The mean intensity has a bell shape similar 
to what is observed for a simulated bar. The se- 
quencing of points also checks that we really mea- 
sure the bar pattern speed. It is worth noting that 
if a linear fit is applied to the plots in Fig. [Ml and 
then the slope is estimated, the pattern speed is 
VL p k, Al km s _1 kpc -1 . 

4.1. Determination of Op by resonances 
analysis. 

We compared ftp obtained from the TW 
method with the one obtained from analyzing 
fundamental resonances of NGC 3367 predicted 
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Fig. 14. — Weighted average position (left) and weighted average velocity (right) along the minor axis. 
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Fig. 15. — Mean intensity along the slits (left) and weighted average velocity as a function of the weighted 
average position (right). 
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by the linear theory. According to the linear the- 
ory at the resonance radii there are expected rings 
of matter with enhanced star formation caused by 
gas shocks. In particular, NGC 3367 has a ring of 
radius ~ 7 kpc formed by spiral arms wrapped by 
more than 180°. Usi ng the diameter of the stellar 
bar 6.7 kpc given by iGarcia-Barreto et al. (2Qo3) 
and assumption that the bar ends near the coro- 
tation, we find that flp = 45 ± 7 km s _1 kpc -1 . 
Here, the errors in fip are due to the rotation 
curve determination and fitting procedure. This 
value agrees within the error limits with the value 
we found by the TW me thod and also accords with 
the v alue reported in Garcia-Barreto Ez Rosadol 



()200 II) of n P = 43 km s -1 kpc" 1 . In Fig. \W\ 



are shown the resonances for NGC 3367 and the 
bar pattern speed determined by the TW method. 
The resonances were obtained from the fit to the 
rotation curve (Fig. [T0|) given in this work. Ac- 
cording to the figure the outer 4:1 resonance is 
located at 7.7 kpc (36.7"), the OLR is at some 
11 kpc (52"), and there is no ILR. However, it 
should be noted that the resulting Qp derived in 
this way should be taken with caution because the 
ring is not necessarily located at corotation. On 
the other hand, given the complex structure of the 
rotation curve in the center and limited resolution, 
it is difficult to judge whether an ILR is absent or 
present. 

5. Discussion 

The main aim of this work is to estimate the 
bar pattern speed of NGC 3367 galaxy using Ha 
kinematic data. For this purpose we first investi- 
gate the sensitivity of results to the data quality on 
numerical galaxy models. We have not intended 
to create an exact numerical model of NGC 3367, 
but rather to use the simulations as toy models. 
For this reason we cannot directly compare the 
lengths and pattern speeds of simulated and ob- 
served bars. We have tested the TW method for 
a simulated galaxy bar against contamination of 
surface density and velocity fields by white noise 
and Gaussian random field perturbations in our 
simulated data. We found that the surface den- 
sity is critical to noise and the bar area should 
have a good signal to noise ratio in order to ob- 
tain reliable results, because it is used as a weight- 
ing function in equation ([1]). In contrast, the ve- 
locity field appears to be less sensitive to errors, 



although they also affect the accuracy of the re- 
sults. This is due to the zero net velocity per- 
turbation introduced in calculus of (V) where the 
errors are partially cancelled out. However, for 
gas-poor, old stellar bars, where the buckling is 
significant, v z component of the velocity field may 
be important. We also checked the influence of the 
posit ion angle on the resulting f2p. In agreement 
with iDebattistal (|2003l ) the results are sensitive to 
PA uncertainties. When the disk major axis and 
the bar become aligned, the errors in Op increase 
faster than in the opposite case. 

Based on the results of the robustness test, we 
apply the TW method to NGC 3367 using pub- 
lished 2D data of Fabry-Perot interferometry. Af- 
ter carefully studying the parameters that could 
affect the results (minimum intensity, range of in- 
tegration, number of strips, high intensity knots 
and arms, position of kinematic major axis), we 
found J7p = 43 ± 6 km s -1 kpc -1 . We also de- 
termine the bar pattern speed by means of reso- 
nances analysis and found ftp = 45 ± 7 km s -1 
kpc -1 . These two results are consistent and in 
agreement with the value report ed previously by 
Garcia-Barreto fc Rosadol (|200lh . The trend in 



variation of fip with the position angle of the kine- 
matical major axis is consistent with the result 
encountered in the simulated barred galaxy. 

As shown in Sect. [2] the variation of the min- 
imum of the surface density for the simulated 
barred galaxy could significantly affect flp deter- 
mination. As we have demonstrated, for NGC 
3367 the adopted minimum of intensity in H a im- 
age is one of the fundamental parameters. Apply- 
ing a clipping allows us to increase the contrast 
of emission in the interarm region, thus improv- 
ing the signal associated with the bar. But, on 
the other hand, if the trend found for the simu- 
lated bar applies for Ha gaseous bar, the clipping 
we used for NGC 3367 could lead to an overesti- 
mation of the fip by more than 50% and the bar 
would lie within a flat region of the curves shown 
in Fig. [8l Thus, if the clipping is not applied, 
our value of Qp indeed can be as small as ~ 15 
km s _1 . Yet, a deeper study using galaxy mod- 
els with complex gas physics is required to verify 
whether the pattern speed is underestimated due 
to clipping. The other two parameters we investi- 
gated essentially restrain the range of integration 
parallel and perpendicular to the kinematic major 
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axis. These are underlying parameters that allow 
us to significantly improve the signal to noise ratio 
near the bar zone. 

When observing the radial velocity field (Fig. , 
an interesting feature can be noted in the nuclear 
region. The isovelocity contours in the center of 
the galaxy that are almost perpendicular to the 
main ones may imply non-circular motions. Since 
the sharp turn in isophotes is frequently associ- 
ated with a secondary bar, the possibility of non- 
circular motions due to different origins is present. 
In particular, the inner gas polar ring or disk, or 
a secondary bar may be responsible for the char- 
acteristic shape of the isovelocity contours. Such 
behavior has been observe d previously in severa l 



galaxies, see for example iMoiseev et al.l ((2004) 



However, due to lack of resolution, in this work 
we were unable to characterize and identify any 
secondary pattern speed. A study with a better 
angular resolution and a comparison with the stel- 
lar counterpart is necessary in order to analyze the 
detailed kinematical features in the center of this 
galaxy. 
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Fig. 16.— Resonances of NGC 3367. The dash- 
dotted horizontal line corresponds to the bar pat- 
tern speed, rip, determined from the TW method. 



6. Conclusions 

In this work we apply the TW method to a sim- 
ulated barred galaxy and test the sensitivity of the 
method to possible sources of error. Then we mea- 
sure the bar pattern speed of NGC 3367 using the 
same method. The results obtained from the ap- 



plication of this method to a simulated galaxy in 
general hold also for the observed one. However, 
the simulations have shown that a significant er- 
ror can emerge due to applying the clipping to the 
surface density. These facts motivated us to as- 
sert some important conclusions within the limits 
of our work. First of all, bearing in mind the prin- 
cipal assumptions of the TW method, we can say 
that it works well enough for the gas phase when- 
ever it is continuously distributed along the bar. 
On the other hand, the validity of applying clip- 
ping to the intensity map should be further inves- 
tigated. We also want to note the importance of 
two dimensional data that allow the exploration 
of a wider range in parameters variation. Given 
the errors of the position angle determination the 
variation of <5PA~ 5° do not affect too much the 
resulting pattern speed for NGC 3367. Finally, 
determining the bar pattern speed by locating the 
corotation resonance, we found that the result is 
similar to that obtained by application of the TW 
method. This fact supports in part the reliabil- 
ity of the method in the case of ionized gas data. 
A deeper study of this galaxy would be worth, in 
particular, the comparison with the pattern speed 
determined from the stellar long-slit spectral ob- 
servations, which would verify the validity of the 
TW method for Ha kinematic data. 
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